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PREFACE 


This  document  reports  the  work  performed  by  McDonnell  Douglas  Astronau¬ 
tics  Company  for  the  Naval  Surface  Weapons  Center,  White  Oak  Laboratory, 
Silver  Spring,  Maryland,  for  the  development  of  an  aerodynamic  code  for 
computing  pressure  loads  on  a  missile  and/or  its  components  at  high  super¬ 
sonic  speeds  and  generalized  attitudes.  The  work  was  accomplished  under 
contract  number  N60921-78-C-0203.  The  Naval  Sea  Systems  Command 
sponsored  the  activity. 

Capabilities  for  analyzing  a  missile  traveling  at  high  supersonic  speeds 
(2.5  <  M^  <  6)  and  generalized  attitudes  have  been  incorporated  into  an 
existing  MDAC  aerodynamic  code  to  provide  the  pressure  loads  at  locations 
specified  by  the  output  of  the  NSWC  computer  codes  PING  and  BING  for 
stress  analysis  by  NASTRAN. 

The  theoretical  background  of  the  method  for  analyzing  the  high  supersonic 
case  is  described  in  this  report.  The  missile  attitude  has  been  generalized 
to  include  the  yaw  and  roll  in  addition  to  the  angle  of  attack.  A  detailed 
description  of  the  input  parameters  is  presented  together  with  test  cases 
and  operating  instructions.  A  complete  listing  of  the  code  is  included. 


CONTENTS 


LIST  OF  SYMBOLS 

vii 

SUMMARY 

ix 

Section  1 

INTRODUCTION 

I 

Section  2 

AERODYNAMIC  THEORY— HIGH 

SUPERSONIC  SPEEDS 

3 

Section  3 

GENERALIZED  MISSILE  ATTITUDE 

3.  1  Coordinate  System  Specifying 

5 

Missile  Attitudes 

3.  2  Effect  of  Yaw  and  Roll  Angles 

5 

on  Aerodynamic  Pressure 

5 

Section  4 

COMPUTER  PROGRAM 

9 

4. 1  Code  Capabilities 

9 

4.  2  Input  Description 

9 

Section  5 

TEST  CASES 

5.  1  Combined  Angle  of  Attack  and 

17 

Yaw 

19 

5.  2  High  Supersonic  Speeds 

21 

5.  3  Test- Case  Input  Cards 

28 

5.4  Test-Case  Output 

5.  4. 1  Combined  Angle  of 

31 

Attack  and  Yaw 

31 

5.4.2  High- Super  sonic  Flow 

50 

Section  6 

CONCLUSIONS  AND  RECOMMENDATIONS 

61 

Section  7 

REFERENCES 

63 

Appendix  A 

PROGRAM  LISTING 

65 

DISTRIBUTION  LIST 

221 

V 


LIST  OF  SYMBOLS 


C 

P 

f 

n 

M« 

X,  y,  z 
Of.  P .  Y 


1  2 

Pressure  coefficient,  (P-In  )  ^  P  U 
Nose  finess  ratio 
Free -stream  Mach  number 
Coordinates 

Angle  defining  missile  attitude  (see  Fig.  1) 


SUMMARY 
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This  report  represents  the  work  performed  under  contract  number  N60921'> 
78-C-0203  for  the  Naval  Surface  Weapons  Center,  White  Oak  Laboratory, 
Silver  Spring,  Maryland,  under  Naval  Sea  Systems  Command  sponsorship 
for  the  development  of  an  aerodynamic  code  for  computing  pressure  loads  on 
a  missile  at  high  supersonic  speeds  auid  generalized  attitudes. 

A  number  of  subroutines  needed  for  computing  the  aerodynamic  pressure  at 
high  supersonic  speeds  (2.5  <  M^<  6)  have  been  developed  and  incorporated 
into  an  existing  MDAC  aerodynamic  code  for  analyzing  the  pressure  loads  on 
a  missile.  For  treating  the  missile  attitude  in  the  general  case,  the  code 
was  extensively  modified  to  reflect  the  required  changes  in  the  input  format, 
code  structure,  and  output  results.  After  the  calculations  are  made,  the 
code  generates  the  pressure  distribution  on  the  missile  at  the  centroid  of 
finite  elements  as  given.  Results  are  produced  as  punched  card  output  in  a 
NASTRAN  bulk-data  format. 


Section  1 
INTRODUCTION 


Development  of  an  aerodynamic  code  (Reference  1)  for  the  calculation  of  die 
pressure  load  on  a  missile  at  the  centroid  of  finite  elements  generated  by 
the  NSWC  computer  code  PING  and  BING  (References  2  and  3)  was  completed 
in  1978  (contract  number  N60921-076-C-0164).  Because  it  was  based  on  die 
numerical  solution  of  a  set  of  linearized  governing  equations  for  die  fluid 
flow,  its  applications  are  restricted  to  the  subsonic  and  low- supersonic 
(M^  <  2.5)  region.  With  recent  improvements  in  missile  technology,  it 
was  deemed  useful  to  extend  its  range  of  validity  to  higher  supersonic  speeds. 
Additionally,  the  limitation  on  missile -attitude  chemges  to  angle -of -attack 
excursions  needs  to  be  removed  to  include  die  attitude  chwges  in  both  yaw 
and  roll.  The  present  effort  was  undertaken  to  meet  these  requirements  in 
order  to  provide  the  necessary  tool  for  structural  analysis. 

To  minimize  the  necessary  effort,  the  earlier  code  was  employed  as  the 
base  into  which  the  high-supersonic  and  the  generalized-missile -attitude 
capabilities  were  incorporated.  Functions  of  the  original  code  such  as 
input  management,  coordinate  specification,  and  transform  and  output 
generation  were  used  with  a  minimum  of  changes. 
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Section  2 

AERODYNAMIC  THEORY  -  HIGH  SUPERSONIC  SPEEDS 


C 


For  high  supersonic  speeds,  the  linearized  method  based  on  small  perturba¬ 
tions  for  lower  Mach  numbers  (M^<  2.5)  can  no  longer  be  employed  for 
analyzing  the  pressure  loads  on  missiles.  Although  exact  numerical  analyses 
such  as  the  method  of  characteristics  or  finite -difference  methods  are  avail¬ 
able,  the  developmental  and  computing  effort  required  for  their  application 
to  the  subject  problem  eliminated  them  from  consideration. 


Because  it  was  considered  to  be  the  most  suitable  for  the  scope  of  the  present 
effort  and  for  possible  future  applications,  the  shock -expans ion  method  was 
selected  for  calctilating  the  pressure  load  on  missiles  at  high  supersonic 
«  speeds.  An  explanation  of  this  method  follows. 

Since  the  Introduction  by  Epstein  (Reference  4)  in  1931,  tiie  shock-expansion 
method  has  been  extensively  developed  by  many  researchers.  The  general¬ 
ized  shock-expansion  method  predicts  the  pressure  distribution  for  flow 
about  pointed  bodies  of  revolution  with  reasonable  accuracy.  In  particular, 
the  shock-expansion  method  (Reference  5)  modified  for  &e  case  where  the 
hypersonic  similarity  parameter  (M^/f^,  ratio  of  free  stream  Mach  number 
to  nose  fineness  ratio)  is  close  to  1,  suits  the  present  application;  its  coded 
version  (Reference  6)  is  used  with  appropriate  modifications. 


The  basic  approach  and  theoretical  derivation  of  the  method  are  described 
in  detail  in  Reference  5  and  will  not  be  repeated  here.  However,  a  brief 
description  of  the  calculational  procedure  is  presented  in  the  following 
sections.  The  method  starts  with  the  cone  solution  for  bodies  of  revolution 
or  the  oblique -shock  solution  for  two-dimensional  flow.  For  pointed  bodies 
of  revolution,  the  pressure  coefficient  after  the  shock  and  for  each 
subsequent  segment  is  obtained  as  &e  sum  of  two  parts.  One  represents  the 
pressure  coefficient  at  zero  incidence  (C  ),  and  the  other  represents 
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the  pressure  coefficient  (C  }  due  to  the  incidence  o .  The  zero-incidence 

pof 

pressure  coefficient  is  computed  using  either  the  second-order  theory  or 
that  developed  in  Reference  7,  depending  on  the  V2d.ue  of  the  parameter 

H,/V 


The  pressure  coefficient  C  for  the  finite  incidence  case  is  calculated  using 

po 

the  following  relations  (References  8  and  9): 


=[AjT  ^  iulJ  - 


*3'"«  lip 


+  [a^t  +A^T/M^^  {ff 


(1) 


where 


A  =  a  +  a,  ,  cos  <|)  +  a,  ,  cos  2<}>  (i  =  1.  2,  ....6) 

X  Uj  X  X^  X  £i p  A 


T  =  sin  Bp  tanp 


The  constants  B  and  a.  ,  are  those  listed  in  Reference  9>  P  and  4>  are  the 

J*  * 

half-cone  angle  and  the  meridian  angle,  respectively. 


The  pressure  coefficient  C  =  C  +  C  obtained  at  each  step  of  computa- 

P  po  p  O’  r  r 

tion  is  used  for  calculating  the  pressure  at  the  end  of  each  segment.  As 
expressed  in  Reference  5: 

P3  =  Pc  -  (Pc  -  (2) 


where 


p-  =  pressure  at  the  end  of  segment. 

^  12 
p^  =  pressure  of  equivalent  cone  =  p^,  C  +Pjp  • 

P2  =  pressure  result  from  either  2-D  expansion  or  compression  between 

segments. 

n  =  function  of  pressure  gradient  and  streamwise  distance. 


For  2-D  flows,  either  oblique  shock  or  expansion  relation  is  used  depending 
on  the  flow -deflection  angle.  The  pressure  load  at  the  centroid  of  each  panel 
is  calculated  by  tracing  the  streamline  over  the  surface  of  the  wing. 
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Section  3 

GENERALIZED  MISSILE  ATTITUDE 


For  specifying  the  attitude  of  a  missile  in  flight,  its  pitch,  yaw  and  roll 
angles  are  usually  required.  An  aerodyneunic  analysis  limited  to  pitch- 
plane  motions  will  be  unrealistic  and  of  rather  restricted  use.  Extending 
the  MDAC  aerodynamic  code  to  include  the  effect  of  yaw  and  roll  expands 
its  range  of  application. 

3.  1  COORDINATE  SYSTEM  SPECIFYING  MISSILE  ATTITUDES 
The  coordinate  system  used  for  defining  the  pitch  angle  o,  the  yaw  angle  P, 
and  the  roll  angle  V  of  the  missile  in  relation  to  the  freestream-flow  vec¬ 
tor  V  is  shown  in  Figure  3-1.  The  successive  rotation  required  to  trans¬ 
form  from  the  flow-field  coordinate  system  (Xj^,  Yj,  Zj^)  to  that  of  the 
missile  (X  ,  Y  ,  Z  )  is  described  as  follows.  First,  the  coordinate  (X., 

Yj^,  Zj^)  is  rotated  about  the  Y  axis  by  an  angle  a  aund  forms  the  intermediate 
coordinates  (X2,  Y^,  Z2).  Next,  it  is  rotated  about  the  Z2  axis  by  an  angle 
P'  to  form  the  coordinates  (X^,  Yy  Z^)*  Finally,  a  rotation  about  (he  X^ 
axis  by  an  angle  Y  leads  to  the  aerodynamic  coordinates  (X  ,  Y  ,  Z  )  for 
aerodynamic  analysis  of  a  missile. 

As  shown  in  Figure  3-1,  either  o,  p  and  Y  or  o,  p'  and  Y  can  be  used  for 
defining  missile  attitude.  For  the  present  case,  the  first  definition  was 
adopted,  i.  e. ,  the  pitch  angle  a,  the  yaw  angle  P  ,  and  the  roll  angle Y  .  Note 
that  both  a  and  P  are  projections  of  angles  in  the  X^-Z,  and  ^j~Y  planes, 
respectively.  In  other  words,  they  are  defined  in  terms  of  angle  in  the  flow 
or  inertial -coordinate  system.  The  reason  for  their  selection  will  became 
clear  when  the  method  for  analyzing  the  pressure  loads  is  discussed. 

3.2  EFFECT  OF  YAW  AND  ROLL  ANGLES  ON  AERODYANMIC  PRESSURE 
When  a  missile  is  considered  as  a  whole  rather  than  as  separate  components, 
accurate  analysis  of  the  aerodynamic  pressure  on  its  surface  in  the  presence 
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Figur«3-1.  Minjie  Attitud*  CoordinatM 


of  yaw  and/or  roll  becomes  somewhat  laborious.  A  simpler  method  that 
provides  results  adequate  for  present  needs  must  be  developed  to  reduce  the 
expected  effort.  For  this  purpose,  a  scheme  consistent  with  the  basic 
approach  used  in  analyzing  the  pressure  load  was  adopted  for  incorporating 
the  additional  effect  of  yaw  and/or  roll. 

As  described  in  Reference  1,  the  approach  employed  is  based  on  the  small- 
perturbation  assumption  which  requires  that  both  angle  of  attack  and  yaw  be 
small.  This  consideration  permits  the  linearization  of  the  governing  fluid- 
dyn2UTiic  equations  and  the  analysis  of  the  pressure  load  on  the  missile  as 
individual  components  (body  and  wings).  Their  interactions,  however,  are 
subsequently  calculated  and  added  to  that  of  the  appropriate  component.  As 
described  below  the  effect  of  yaw  and  roll  on  the  aerodynamic  load  on  the 
missile  is  similarly  analyzed  taking  advantage  of  this  approach. 
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The  combined  effect  of  angle  of  attack  a  and  yaw  p  on  the  body  can  be 
viewed  as  that  of  an  equivalent  angle  of  attack  o'  with  a  rotation  of  the  body 
about  its  axis  by  an  angle  Y*.  From  the  geometry,  one  has 


Ofi  =  arc  cos  [l.  O/  V  1  +  tan^  a  +  tan^p  j 

(3) 

Y'  =  arc  tan  j^tanp/tanoj 

(4) 

As  defined,  the  effect  of  roll  V  represents  only  a  shift  in  the  angular  distri¬ 
bution  of  the  pressure  load  on  the  body  and  is  incorporated  by  substituting 
the  algebraic  s\am  of  Y  and  Y'  or  Y"  for  Y'  in  the  analysis. 

The  combined  effect  of  angle  of  attack  o,  yaw  p  and  the  roll  Y  on  the  wing 
pressure  load  is  incorporated  by  an  apparent  angle  of  attack  o  '’  which  is 

O'"  =  arc  sin  ^sinY"  sin  p'  cos  a  +  cos  Y"  sin  oj  (5) 

where 


The  advantage  of  the  present  approach  for  incorporating  the  effect  of  yaw  and 
roll  is  that  it  allows  maximum  use  of  the  output  of  the  existing  analysis  and 
consequently  minimizes  the  effort  required  to  develop  a  more  versatile  code. 
The  practical  necessity  of  code  reorganization,  core-storage  rearrangement, 
and  appropriate  assembly  of  additional  subroutines  have,  however,  been 
time  consuming,  and  repeated  computer  runs  were  required  for  program 
checkout. 
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Section  4 

COMPUTER  PROGRAM 


The  MDAC  aerodynamic  code  (Reference  1)  previously  developed  for  NSWC 
under  the  contract  number  N60921-76-C-0164  has  been  extensively  modified 
to  incorporate  the  high-supersonic -speed  and  the  generalized  missile- 
attitude  capabilities.  Using  the  output  of  NSWC  computer  codes  PING  and 
BING  which  specify  the  finite -element  and  grid -point  distribution  on  the  mis¬ 
sile,  the  present  code  calculates  the  desired  pressure  load  at  the  centroid  of 
each  element.  At  the  user's  option,  the  pressure  load  may  be  generated  as 
punched-card  output  in  the  bulk-data  format  of  NASTRAN  for  subsequent 
structural  analysis. 


4.  1  CODE  CAPABILITIES 

The  present  code  as  modified  is  capable  of  calculating  the  pressure  load  for 
a  missile  or  its  components  (body  or  wing)  at  either  subsonic  or  supersonic 
speeds.  The  pitch -plane  restriction  on  its  attitude  has  been  removed  to  allow 
both  yaw  and  roll  in  addition  to  pitch-plane  motions.  Due  to  the  small- 
disturbance  assumption,  the  code  is  valid  for  slender  bodies  of  circular  or 
near  circular  cross-section  at  small  angle  of  attack  or  yaw. 

The  code  requires  a  field  length  of  150,  000  when  used  on  the  CDC  CYBER  174 
computer  with  KRONOS  Z,  1  system.  For  subsonic  or  low-supersonic  cases, 
approximately  3  minutes  computing  time  is  required.  For  higher  supersonic 
(M  >  2.5)  cases,  much  less  computing  time  («  1.  0  min. )  is  required. 


4.  2  INPUT  DESCRIPTION 


Input  parameters  for  numerical  computation  have  been  modified  as  needed, 
based  on  the  MDAC  aerodynamic  code  previously  developed.  Only  minor 
changes  were  required  for  specifying  either  the  high-supersonic-speed  case 
or  the  generalized -attitude  case.  However,  for  user's  convenience,  a  com¬ 
plete  list  of  input  parameters  has  been  included  in  the  following  description: 
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•~PitE€EBim  MOB  NOT  FUMED 
BLANK 


XMACH 

PINF 
DA  DEG 

SID 


YAW 

BRCJLL 


Mach  number  of  flow.  Higher  supersonic  routine 
employed  when  greater  than  2.  5,  default  value  =  2.  0. 
Free  stream  pressure,  default  value  =  14.7  psi. 
Incremental  angle  of  attack  of  missile  in  degrees, 
default  value  =  0. 

Case  identification  number  for  each  missile  com¬ 
ponent  (up  to  20  allowed);  default  values  have  been 
set  at  1.  Also  used  for  selecting  the  sign  conven¬ 
tion  for  pressure  output.  Code  uses  the  accepted 
rule  for  pressure,  i.  e.,  positive  acting  toward  the 
surface  and  negative  acting  away  from  the  surface. 
For  opposite-sign  convention,  set  SID(L)  as  nega¬ 
tive  for  pressure  output  on  the  desired  component 
L  (in  the  order  of  body,  wing  surfaces,  upper  or 
lower). 

Angle  of  yaw  in  degrees  with  respect  to  the  flow. 
Missile  roll  about  body  axis  (degree). 


The  following  input  parameters  are  required  for  the  case  where  missile 
body  is  included: 


» 


NBPFL 

XB  <  51 

RB  <  51 

ZD  <  51 
ARB 


Number  of  body  profile  stations  for  specifying  body 
geometry  inputs  r  =  r(x);  up  to  51  stations  allowed. 
Axial  coordinate  of  body  stations  in  the  direction  of 
flow  for  specifying  the  body  profile. 

Radial  coordinate  of  body  profile  at  corresponding 
body  stations. 

Body  camber  at  corresponding  body  stations. 

Angle  of  attack  of  missile  body  in  degrees,  default 
value  =  0. 


As  is  similar  to  the  overall  organization  of  the  code,  the  input  parameters 
and  data  have  been  organized  into  two  main  groups  in  accordance  with  their 
assigned  functions.  The  first  group  contains  all  the  parameters  needed  for 
the  specification  of  case  selection,  missile  configuration,  and  flow  conditions. 
They  form  the  input  parameter  of  the  namelist  UNIFID.  The  second  group 
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specifies  fhe  finite -element  and  grid-point  system  as  generated  by  the  NSWC 
computer  codes  BING  and/or  PING.  It  includes  the  specification  for  the 
coordinate  system  used  for  the  finite  elements. 

For  the  convenience  of  the  user,  card  format  for  inputs  are  described  in 
their  sequential  order  as  required  by  the  code. 


A.  Aerodynamic  Control  Parameters. 

Card  1  Title  and  case  description,  format  (20A4),  alphanvimeric. 
Card  2  $UNIFID,  begins  at  Column  2  with  the  symbol  $. 

Card  3  And  additional  cards  if  necessary,  contain  the  following 
as  input: 


ICASE  =  1 
=  2 
=  3 

NTRANS 


IPUNCH  =  0 
=  1 

IRW  =  0 

=  1 

=  2 


POLAR 


Wing  only. 

Body  only. 

Complete  missile. 

Number  of  coordinate  transformations  required  to 
convert  the  coordinates  for  the  finite  element  to 
the  coordinates  used  in  the  aerodynamic  code.  Two 
transformations  are  allowed  for  each  missile  com¬ 
ponent,  in  the  following  order;  body,  wings,  and 
nacelles. 

No  punched  card  output,  default  value. 

Punched  card  output  requested. 

Computation  results  are  not  saved  on  tape,  default 
value. 

Aerodynamic  computation  results  are  saved  on 
tape  12  for  restart,  first  run  only. 

Restart  run,  bypass  all  aerodyncunic  computations 
and  begins  computation  at  the  start  of  pressure 
interpolation. 

Number  of  incremental  angles  of  attack,  default 
value  =  0. 


The  following  inputs  are  required  for  the  case  where  the  missile  wing  is 
included: 
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ISOLID  =  0 


Wing  of  built-up  construction,  pressure  load  on 
both  upper  and  lower  wing  surface  generated, 
default  value. 


=  1 

Wing  of  solid  section,  sum  of  the  pressure  loads 
on  wing  surface  generated. 

IFORM  =  0 

Wing  of  built-up  construction,  both  upper  and  lower 
surface  height  (ZWI)  are  given  at  identical  locations 
(XWI,  YWI). 

=  1 

Wing  of  built-up  construction,  upper  and  lower 
surface  height  (ZWI)  are  not  given  at  identical 

locations. 

=  2 

Flat  wing  surface. 

NWPI 

Number  of  coordinate  points  specify  the  wing  sur¬ 
face  contour  (<30).  Code  allows  a  maximum  of 

20  wing  surfaces  (upper  and  lower).  Each  wing 

surface  must  be  defined  in  the  form  of  equation  of 

a  surface,  ZWI  =  f(XWI,  YWI),  A  maximum  of 

30  points  are  permitted  with  the  provision  that  die 

first  four  points  are  restricted  to  the  specification 
of  the  corners  of  the  wing  only.  The  coordinates 
for  the  corners  must  be  given  in  the  following  order. 
Starting  with  the  innermost  point  of  the  leading 
edge,  the  remaining  corners  are  to  be  given  in  a 
clockwise  order  viewing  from  the  above. 

NWING 

Number  of  wings  or  tails.  A  meucimum  of  10  is 
allowed,  default  value  =  0. 

XWI(i,  j). 

Wing-surface  coordinates  where  indices  i  and  j 

YWI{i,  j). 

designate  the  surface  coordinate  point  and  wing 

ZWI(i,  j) 

surface,  respectively. 

DIKED 

Dihedral  angle,  degrees. 

PIVOT 

Indicates  dihedral,  0  for  no  dihedral,  and  >0  for 

• 

dihedral  present. 

ARW 

Angle  of  attack  of  wing  in  degrees,  with  respect  to 

—  1 

the  body  aucis  if  complete  missile  is  considered, 

default  value  =  0. 

rV 
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■  i 

The  following  input  parameters  are  required  for  the  complete -missile  case 
(ICASE  =  3)  except  for  high-supersonic-speed  (M  >  2.  5)  case: 


NWPANL 

XLG 

XTG 


Number  of  wings  directly  connected  to  the  body. 
Axial  distance  at  the  intersection  of  wing  leading 
edge  and  body  for  each  wing. 

Axial  distance  at  the  intersection  of  wing  trailing 
edge  and  body  for  each  wing. 


Card  4  $END 


Begin  at  column  2  the  symbol  $  for  ending  input. 


B.  Finite -Element  and  Grid  Inputs,  The  second  group  of  input  cards 
are  arranged  immediately  following  Card  4.  as  follows: 


Cards  5  and  6  Coordinate -system  specifications  in  the  format  of 

NASTRAN  bulk-data  deck  [page  2.4-49  -  2.4-54, 
NASTRAN  User's  Manual,  NAS  SP-22 1(0 1),  1972]. 
Three  position  vectors  —  A,  B,  and  C  —  are  used 
to  define  the  coordinate  system.  The  first  defines 
the  origin,  the  second  defines  the  Z  axis,  and  the 
third  defines  a  point  in  the  XZ  plane. 


For  Card  5  the  following  format  is  used: 

Col  1-8,  Coordinate  system  CQJRDJRbb  for  rec¬ 
tangular,  C0RDJCbb  for  cylindricad,  and 
CQlRDJSbb  for  spherical,  where  J  identi¬ 
fies  the  coordinate  system  numbers. 

Col  9-16  Coordinate  identification  number  J,  integer. 
Col  17-24,  Reference  coordinate  system,  integer, 
optional. 

Col  25-32,  33-40,  41-48,  Components  of  vector  A 

(3F8.2). 

Col  49-56,  57-64,  65-72,  Components  of  vector  B 

{3F8.2). 


For  Card  6: 


Col  1-8,  blank. 

Col  9-16,  17-24,  25-32,  Components  of  vector  C 

(3F8.2). 


iMeoowMtfi.4.  aauat.Mi 
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Cards  5  and  6  are  to  be  repeated  as  many  times  as  there  are  coordinate 
systems  (NTRANS)  required  for  specifying  the  finite  element.  They  are  to 
be  arranged  in  the  order  of  body  and  wings. 


The  finite  element  and  grid  specifications  are  arranged  immediately  follow¬ 
ing  the  coordinate -system  input  in  the  same  order,  i.  e. ,  body,  wings. 

Card  7:  Specifies  the  number  of  grid  cards  (NGRIDP)  and 

element  cards  (NE),  Format  (ZllO). 

The  grid  and  the  element  cards  as  generated  by  die  NSWC  computer -code 
BING  are  sorted  into  two  separate  groups.  They  are  then  placed  immedi¬ 
ately  after  card  7  with  the  grid  cards  in  front. 

Some  remarks  regarding  the  input  preparation  need  to  be  made  here  to 
facilitate  the  use  of  the  aero  code. 

For  the  body-only  case,  inputs  of  the  body  profiles,  expressed  as  r  «  r(x), 
for  aero-panel  generation  must  be  in  the  aero -coordinate  system.  It  is 
required  that  x  be  the  axial  distance  from  the  tip  and  r  the  radius  of  the  body 
cross  section  (Figure  4-1).  For  the  finite  element  and  grid  point,  their 
relationship  to  the  basic  coordinate  system  and  the  basic  coordinate  system 
to  aero -coordinate  system  are  to  be  specified  in  the  coordinate  input  card 
immediately  following  the  namelist  input.  The  grid  card  and  the  element 
card  are  to  be  separated  and  placed  after  the  coordinate  system  specification, 
with  a  card  specifying  their  numbers  placed  behind  die  coordinate  card. 

The  input  for  the  wing-only  case  requires  the  specification  of  both  the  wing 
planform  and  its  surface  contour.  The  input  for  the  description  of  the  wing 
planform  requires  the  coordinates  of  corner  points  and  break  points.  For  each 
individual  wing,  inputs  for  those  of  quadrilateral  shape  are  their  four  corner 
coordinates.  The  codes  have  allocated  the  first  four  values  of  the  coordinates 
XWI,  YWI  and  ZWI  for  each  wing  or  wing  region,  arranged  in  a  counterclock¬ 
wise  order  starting  from  the  leading-edge  corner  of  break  point.  In  the  case 
of  triangular  wings,  there  are  only  three  corner  points.  To  satisfy  the  same  ■ 
requirement,  the  apex  or  wing  tip  is  to  be  considered  as  two  coincident 
points  of  the  same  coordinate.  The  inputs  are  to  be  arranged  like  those  of 
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Figurt4-1.  Body  Profiia  Input 

the  quadrilateral  shape,  with  leading-edge  body-corner  point  as  the  starting 
point.  Wings  with  complex  shapes  (see  Figure  4-2)  must  be  subdivided  into 
appropriate  quadrilateral  or  triangular  regions  and  the  inputs  prepared 
accordingly.  The  only  requirement  is  that  the  control  chords,  apart  from 
either  the  leading  edge  or  the  trailing  edge,  must  be  in  the  streamwise 
direction. 


All  inputs  for  the  wing  surface  are  to  be  given  in  the  form  of  the  equation  of 
surface  Z  =  f(x,  y)  in  terms  of  XWI,  YWI  and  ZWI,  Starting  with  the  fifth 
value,  XWI  (5, 1),  YWI(5, 1),  ZWI(5, 1),  their  total  number  may  not  exceed 
26  points.  For  wings  with  flat  surfaces,  only  the  first  four  points  are 
necessary. 

Inputs  defining  wing  thickness,  camber  and  twist  have  been  eliminated  to 
reduce  the  input  preparation  effort.  The  code  generates  the  required  values 
of  camber,  thickness  slope,  and  twist  angle  for  the  wing  using  the  existing 
surface -fit  scheme  from  the  input. 

Similar  to  the  body-only  case,  the  finite -element  and  grid -point  coordinate - 
system  cards  are  placed  immediately  behind  the  namelist  input.  They  are 
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YWI(I^) 

YWKM) 

ZWI(1,3> 

k 

XWI(2.3I  XWI(34i,  XWI(44I 
YWt(24)  YWII34>.  VWII44I 
ZMI(24)  ZWI(34).  ZWII44I 


NWING  >  3,  NTRANSdl  =  1,  1.  1, 

If  Built-Up  Wing  it  Used,  NWPKl),  NWPI(2),  and  NWPIOI  Mutt  be  Specified 


Figura4-2.  Subdivision  of  Complax  Wing 


« 


followed  by  a  card  specifying  the  number  of  grid  cards  and  element  cards. 

The  grid -point  and  element  card  as  generated  by  PING  occupies  the  last  posi¬ 
tion  in  the  input. 


The  coordinate  system  used  for  the  wing  can  be  any  one  of  three  systems, 
i.e.,  rectilinear  (x,  y,  z),  cylindrical  (r,e,z)  and  spherical  (r.G  In  terms 

of  the  input  parameters,  they  are  represented  by  XWI,  YWI  and  ZWI,  respec¬ 
tively,  for  all  three  systems. 

For  the  complete -missile  case,  the  input  can  be  considered  as  a  combination 
of  these  two  cases,  i.e.,  the  body-only  and  the  wing-only  cases.  Only  the 
parameters  ICASE  and  NTRANS  need  be  corrected  to  reflect  the  changes. 

Both  coordinate -system  cards  and  finite -element  cards  are  to  be  grouped 
together  in  the  order  of  body,  wings,  and  nacelles. 

As  a  result  of  the  approach  used  for  the  analysis,  the  computations  for  the 
high-supersonic -flow  case  are  considerably  simpler.  Accordingly,  the  code 
treats  the  complete  missile  as  a  rule. 


A9COOfVMl.l.  OOC/OA.A] 
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Section  5 
TEST  CASES 

Test  cases  were  selected  for  aerodynamic  analysis  by  the  present  code  for 
the  purpose  of  verifying  the  theoretical  approach  and  code  modifications.  In 
Reference  10,  the  experimental  results  from  tests  conducted  for  a  model 
with  delta  planform  wings  and  cylindrical  body  at  supersonic  speeds(M^=  2.  0 
and  2.  58)  are  available  and  suitable  for  our  purposes.  The  model  configure • 
tion  with  dimensions  is  shown  in  Figure  5-1.  Location  of  the  pressure  taps 
for  wing  and  body  are  shown  in  Figures  5-2  and  5-3,  respectively.  Tests 
were  made  for  angles  of  attack  and  yaw  varying  from  0“  to  5”. 


t0.7t8  in. 


3.403  In. 


6.803  in. 


4.137  bt. 


Figura  5-1.  Taft  Modal  (Rafaranea  10) 
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To  demonstrate  the  ability  of  the  code  in  analyzing  the  pressure  load  on  a 
missile  at  high-supersonic  speed  and  generalized  attitude,  computations  were 
made  for  missiles  flown  at  (a)  combined  angle  of  attack  (5°)  and  yaw  (5**)  at 
Mach  2.0  and,  (b)  high-supersonic  case  at  Mach  2.  58  and  5°  angle  of  attack. 


5.  1  COMBINED  ANGLE  OF  ATTACK  AND  YAW 
The  following  is  a  list  of  input  cards,  in  sequence: 

(I)  Title  Card. 

Descriptive  title  of  the  case. 

FORMAT:  (20A4). 

(II)  Namelist  UNIFID  input. 

Aerodynamic  control  parameters  and  geometries  of  missile. 
FORMAT:  Real  or  integer  in  accordance  to  the  variable. 


Test  Case: 


ICASE  =  3 
NTRAN(1)=1,  1 
IS0LID  =  1 
NBPFL  =  15 
NWPI(1)=4 
NWINGsl 
IF0RM(1)=2 
ARB  »  5.0 


Complete  missile. 

No  intermediate  coordinate  system  for  wing  or  body. 
Solid  wing,  sum  of  upper  and  lower  pressure. 
Number  of  input  points  for  body  profile. 

Number  of  input  points  for  wing  surface. 

One  wing  section. 

Flat-wing  surface. 

Angle  of  attack  of  missile. 


YAW  =  5.0 
XMACH  =  2.0 
NWPANL  =  1 
XLG(1)=6.053 
XTG(1)=10.  718 


Angle  of  yaw  of  missile. 

Flow  Mach  number. 

Number  of  wings  directly  connected  to  the  body. 
Axial  distance  of  the  inner  leading  edge. 

Axial  distance  of  the  iimer  trailing  edge. 


For  body  profile ; 

XB(1)=0. 0,  0.5,  1.0,  1.5,  2.0,.... 

RB(1)=0. 0,  0.0881,  0.  1763,  0.2664,  0.3526,.... 


For  the  wing  surface,  only  four  points  are  needed  for  a  solid  wing: 
XWI(1,  1)=6.  503,  10.718,  10.718,  10.718. 

YWI(1,  1)=0.6,  0.6,  4.  137,  4.  137. 

ZWI(1,  1)=0.  0,  0. 184,  0.  0,  0. 184. 
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(lU)  Coordinate  Card. 

Notes:  The  coordinate  cards  are  needed  for  input  of  vectors  A,  B,  and 

C  which  are  used  in  defining  the  relation  between  the  coordinate 
systems  of  the  finite  element  (X,  Y,  Z)  and  the  aero  codes  (X^, 

Y.,1  Z-).  In  accordance  with  NASTRAN,  vector  A  defines  the 

CL  <* 

origin,  vector  B  defines  the  Z  axis,  and  vector  C  defines  the 
X-Z  plane. 

FORMAT: 

NASTRAN  bulk-data  deck,  first  card  (A8,  218,  6F8.2)  and 
second  card  (A8,  3F8. 2). 

Test  Case: 

Two  sets  of  coordinate  cards  are  required,  one  for  the  body 
and  one  for  the  wing.  The  body  finite -element  coordinate  sys¬ 
tem  has  its  origin  coincident  with  the  aero  code  and  its  Z  axis 
coincident  with  the  X  axis  of  the  aero  code  as  shown  in  the 
sketch.  Hence,  the  vector  A  has  component  (0.0,  0.0,  0.0), 
vector  B  has  component)(l.  0,  0,  0,  0.  0),  and  vector  C  has 
component  (1.  0,  0.0,  1.0). 

The  finite -element  coordinate  system  for  the  wing,  as  shown 
in  the  sketch,  requires  the  input  of  vector  A  with  component 
(0.  0,  0.  0,  0.  0),  vector  B  (0.  0,  0.  0,  1.  0),  and  vector  C 
(1.0,  0.0,  1,0). 

(IV)  Finite -Element  Cards. 

Notes;  The  finite -element  and  grid-point  cards  as  generated  by  BING 

and  PING  are  used  for  defining  element  coordinates.  They  must 
be  divided  into  two  separate  groups  with  the  grid-point  card 
preceding  the  element  card.  A  card  specifies  that  the  numbers 
of  grid  card  and  element  card  in  a  format  (2110)  must  be  placed 
before  them. 

FORMAT; 

For  grid  and  element  number  card  (21 10),  grid  and  finite- 
element  cards  must  be  the  same  as  generated  by  BING  or  PING. 
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Test  Case: 


As  shown  in  the  input  listing,  first  set  of  finite -element  grid 
points  and  element  cards  are  that  of  the  body  and  the  subse¬ 
quent  set  is  for  the  wing  surface. 
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Results  of  Computations 

Using  the  inputs  as  described,  computations  were  made  for  the  pressure  load 
on  the  body  and  wing  of  the  missile  tested  in  Reference  10.  In  Figure  5-4, 
the  pressure  load  on  the  body  at  the  meridian  angles  of  45°  and  135°  were 
compared  with  that  of  the  experiment.  Good  correlation  with  the  experimental 
measurements  was  obtained.  In  Figure  5-5,  comparisons  with  the  experi¬ 
ment  were  made  for  the  pressure  load  on  the  wing  and  the  results  were 
satisfactory. 

It  is  seen  that  the  code  can  be  used  for  analyzing  the  aerodynamic  loads  on 
missile  at  attitudes  that  may  include  yaw  and  roll.  Its  applicability  is  there¬ 
fore  considerably  extended. 

5.2  HIGH  SUPERSONIC  SPEEDS 

The  input  for  this  case  differs  fr(»n  the  previous  case  only  in  the  Mach 
number,  i.  e.,  2.58  instead  of  2.0.  All  the  remaining  input  parauneters  have 
been  retained  for  this  computation  and  hence  will  not  be  repeated  here. 
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Fig.  54.  Pressure  Distribution  on  Body  at  Mach  2.0, 5°  Angle  of  Attack  and  5^  Yaw 
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Fifura  5-5.  Pranur*  DittrMwtion  on  Wing  at  Mach  2.0, 9*  Angla  of  Attack  and  5**  Yaw 
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Computational  Results 

Results  of  computations  for  flow  at  Mach  2.  58,  illustrated  in  Figure  5-6  for 
the  body  and  Figure  5-7  for  the  wing,  show  that  the  pressure  distribution  on 
the  body  correlates  poorly  with  the  test  data.  It  is  believed  that  this  results 
from  a  substantial  effect  of  the  wing  on  the  body.  This  effect  should  diminish 
with  increasing  Mach  number.  The  pressure  distribution  on  the  wing  surface 
appears  to  be  acceptable  over  most  of  the  wing  surface. 

It  should  be  pointed  out  that  due  to  the  scarcity  of  experimental  data  for 
pressure  distributions  on  complete  missiles,  the  test  case  described  was 
computed  for  a  Mach  number  value  slightly  above  the  lower  limit  for  the  high- 
supersonic -speed  analysis.  However,  for  bodies  of  simple  geometries 
(circular  cones  and  triangular  wing),  the  test  cases  were  computed  at  Mach 
numbers  6.  85  and  5.  05,  and  the  results  correlated  rather  convincingly.  In 
Figure  5-8,  pressure  coefficients  on  circular  cone  with  half-cone  angle  of 
12.  5°  at  angles  of  attack  of  3°,  6“,  and  12“  were  calculated,  and  the  results 
were  very  good  when  compared  with  available  data.  In  Figures  5-9  and  5-10, 
pressure  coefficients  and  shock  shape  for  an  ogive  at  an  e  of  attack  were 


Expcrimtnt  (Rafcrenct  10) 


Figure  SS.  Preuure  Oiitribution  on  Body  at  Mich  2.58, 6®  Angle  of>Kttack 
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Figure  5-7,  Presaure  Distribution  on  Wing  at  Mach  2.58, 5**  Angle  at  Attack 

computed,  and  the  results  were  satisfactory  as  shown.  In  Figure  5-11,  the 
pressure  on  a  75®  sweep,  triangular  wing  was  calculated  using  the  shock- 
expansion  method.  Comparisons  with  data  show  good  correlation. 


o<>uai.A 
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M_  -  6.86 

Half  Cone  Angla^  •  12.6° 


0.4 


0“  30*»  80"  90P  120“  IMP 

MarMlan  Angla# 


Figura  54.  Pranurt  Coafficiant  for  Cona 
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Zero  An0l«of 
Attack 


O  Exparimant  (Rafaranca  5) 
~  Computed 


5.  3  TEST -CASE  INPUT  CARDS 


AfCOOATIVMi.  OOC/OI^l 


28 


tiriKIIIENT  CASE  rci  CNECRINE  lENEtALIlEA  MISSILE  ATTITUSE  M«»E 
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5.  4  TEST-CASE  OUTPUT 

5.  4. 1  Combined  Angle  of  Attack  and  Yaw 
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Section  6 


CONCLUSIONS  AND  RECOMMENDATIONS 

The  MDAC  aerodynaunic  code  has  been  extended  to  treat  (a)  high  supersonic 

speeds  (M  >  2.5)  and  (b)  generalized  missile  attitudes.  As  a  result  it  is 

00 

capable  of  analyzing  the  missile  under  more  realistic  flight  conditions. 

Effort  for  input  preparation  remains  the  same  for  all  practical  purposes 
because  only  a  few  inputs  were  added  for  the  specification  of  attitudes.  Com¬ 
puting  time  for  high  supersonic  flow  analysis  will  be  only  a  small  fraction  of 
that  for  the  low-speed  case.  The  original  flexibility  in  the  coordinate  system 
has  been  maintained  for  the  convenience  of  users.  Outputs  are  produced  as 
computer  printout  and  punch  cards  (if  desired)  in  the  format  of  NASTRAN 
bulk  data. 

With  these  additional  capabilities  incorporated,  the  code  represents  a  versa¬ 
tile  tool  for  aerodynamic  analysis  of  missiles  and  forms  a  unique  link  between 
the  flight  dynamics  and  structural  analysis  of  the  missile  under  a  wide  variety 
of  flight  conditions  that  it  may  encovinter  in  actual  operation.  Furthermore, 
the  large  effort  necessary  for  preparing  the  input  to  the  NASTRAN  is  reduced 
when  it  is  used  in  combination  with  the  NSWC  computer  codes  PING  or  BING. 

For  missile  flight  at  high  supersonic  speeds  (2  <  8)  the  problem  of  aero¬ 

dynamic  heating  becomes  important  and  needs  to  be  analyzed.  As  an  example, 
for  an  insulated  flat  plate  aligned  parallel  to  the  flow,  the  surface  tempera¬ 
ture  can  reach  anywhere  from  1000**R  to  7000°R  for  a  speed  range  of 
M^=  2  to  8.  Additionally,  during  maneuvering,  severe  heating  may  occur  as 
a  result  of  shock  impingement  on  the  missile,  e.  g. ,  nose  bow  shock-wing 
interaction.  Thus,  the  induced  thermal  stress  can  seriously  affect  the  strucv 
tural  integrity  of  the  missile,  and  failure  of  a  missile  component  can  result. 
Consequently,  development  of  capability  for  analyzing  the  aerodynamic  heat¬ 
ing  to  complement  the  pressure-load  calculation  will  be  an  appropriate  task 
and  is  recommended  as  a  future  effort.  The  required  flow-field  quantities 
such  as  the  pressure,  velocity,  and  density  can  be  derived  conveniently  from 
the  present  code  as  constructed.  This  feature  can  be  expected  to  minimize 
the  necessary  work  compared  to  that  of  an  independent  effort. 
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